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Performance of Graphene Thermal Photon Detectors 
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We analyze the performance of graphene microstructurcs as thermal photon detectors and deduce the range of 
parameters that define a linear response. The saturation effects of a graphene thermal detector that operates 
beyond the linear range are described in detail for a single-photon detector (calorimeter). We compute 
the effect of operating beyond this linear range and find that sensitive detection occurs for such non-linear 
operation. We identify the optimum conditions and find that single-photon detection at terahertz (THz) 
frequencies should be feasible. 



I. INTRODUCTION 



Modern photon detectors are widely employed 
in sensitive applications ranging from astrophysical 
observations^'^ to quantum communications."^ Some of 
the most sensitive detectors employ a thermal detec- 
tion mechanism, where the photon energy (or power) 
is converted to a temperature increase. This tempera- 
ture increase is sensed through a change of the electri- 
cal properties, such as a resistance change for a super- 
conducting device biased on its superconducting transi- 
tion. Transition-edge sensors (TESs)^ are well developed 
and their sensitivity can approach the fundamental lim- 
its for thermal detectors.^ Hot-electron TES detectors, in 
which only the electrons are heated, have achieved single- 
photon sensitivity in the near-infrared^ and the mid- 
IR.^'* Other detectors (such as quantum capacitance^ 
and kinetic inductance^" detectors) are also being de- 
veloped for these applications. 

Lower-dimension carbon nanostructures (graphene 
and individual carbon nanotubes) are being considered as 
thermal photon detectors. ^^"^^ Since these systems have 
very few electrons that are active in the conduction pro- 
cess, the photon energy can cause significant heating. In 
addition, the very stiff carbon lattice should ensure that 
the clcetron-phonon interaction is weak, so cooling by 
phonon emission is small. 

In this paper we investigate the limits on the sensi- 
tivity of graphene-based thermal detectors which employ 
a Johnson-noise temperature readout. Neither the dy- 
namic range of detector operation nor the effect of oper- 
ating beyond the linear range of thermal response have 
been studied previously. We concentrate on graphene, as 
carbon nanotubes have large quantum contact resistance 
and small shunting capacitance of the contact, which 
makes THz photon coupling a challenge. Also, graphene 
has been shown to have simple Drude behavior of the 
conductivity from DC into the THz range, ^® so modeling 
the absorption is feasible. The graphene detector exhibits 
saturation of the response when operated beyond the lin- 
ear range, as defined below. We conclude that, with a 
Johnson-noise readout, operation beyond the linear range 
(near-equilibrium) provides the most sensitive detection 
and a graphene calorimeter can detect individual THz 
photons. Our study shows the limits and constraints of 
the Johnson- noise temperature readout. In this paper we 
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FIG. 1. Schematic of bolometer or calorimeter. C is the heat 
capacity and G is the thermal conductance. In a graphene de- 
tector, C and the thermometer function are provided by the 
electron subsystem. G is provided by the coupling to the con- 
tacts and substrate as well as emitted Johnson noise, (b) Tem- 
perature response, T{t), for a calorimeter. For graphene, re- 
sistance is nearly independent of temperature, (c) Schematic 
of antenna coupling to a small graphene device at the center. 
For a THz detector, the antenna linear dimension is a few 100 
/im. The low frequency contacts are not shown. 



first introduce the qualities of a desirable calorimeter and 
summarize the thermal properties of graphene. Then we 
describe how accurately one can read the temperature of 
graphene, both when the system is near equilibrium and 
after absorption of a 1-THz photon. Finally, we assess 
the use of graphene THz photon counter. 



II. GRAPHENE AS A CALORIMETER 

A calorimeter works by absorbing a photon and as a 
result, changing its temperature, which is then read out 
to determine the photon energy. Such a thermal detector 
can also be used to measure the incident power due to 
a flux of photons, in which case it is called a bolometer. 
Whether it is optimal for power detection to operate in 
the integrating mode or in the photon counting mode is 
determined by the detector response time and the pho- 
ton arrival rate, provided the detector can detect a single 
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photon.^ The figure of merit for a calorimeter is the en- 
ergy resolution of the detector - the ratio of photon en- 
ergy to the fluctuations in energy during a measurement. 
In the linear range of response for a thermal detector, 
the intrinsic resolution is not affected by the size of the 
input signal. We define this linear range by requiring 
that the initial temperature increase due to a single ab- 
sorbed photon (Fig. lb) be AT = (Tpcak-To) < Tq; Tq 
is the device temperature with no photons and Tpeak is 
the electron temperature shortly after photon absorption. 
For devices based on graphene, we require AT < aTo to 
specify the linear range, and in this work we set the limit 
a 0.1 as a reasonable choice. 

Before discussing the performance of graphene as a 
photon detector, we first review the thermal properties 
of graphene. At low temperatures the heat capacity, 
C, is due to the electron specific heat and scales with 
temperature as C = 'yAT, where 7 is the Sommerfeld 
coefficient of the electron system, and the thermal con- 
ductance scales as Gcph « AT^ for phonon cooling; A 
is the graphene area.^^ These results are for finite elec- 
tron densities away from the Dirac point in the stan- 
dard model. Cooling by low-frequency photon emis- 
sion can also contribute to the thermal conductance as 
Gphoton ~ k^B, where B <^ ksT/h is the coupled band- 
width to an impedance-matched load.^^ In our case, we 
consider graphene which is impedance matched to an am- 
plifier. Thus, there is no amplifier backaction (no electro- 
thermal feedback). We assume that electron-electron 
interactions rapidly convert all the photon energy E to 
electron thermal excitations at an initial temperature 
Tpoak above Tq. Indeed, very fast electron-electron en- 
ergy exchange for graphene has been reported. 

If outdiffusion of heat through the contacts^^ is sup- 
pressed, the thermal conductance G is due to the emis- 
sion of phonons and microwave photons^^ by the heated 
electrons: 

G = Gcph + Gphoton- (1) 

Suppression of diffusion cooling out the contacts is de- 
sirable for a high sensitivity detector^^'^^ to allow a 
longer averaging time of the output signal. This can 
be achieved by transparent contact between graphene 
and superconducting leads, ^^"^^ where the superconduc- 
tor energy gap prevents the outfiow of thermal excita- 
tions from the graphene.^'' Opaque (tunneling) super- 
conducting contacts (superconductor-tunnel insulator- 
graphene) also prevent heat outdiffusion. ^^^^^ With such 
contacts, the THz coupling can be efficient due to the 
finite tunnel junction capacitance. Additionally, the low- 
frequency resistance is temperature dependent, which 
provides a resistance readout of the graphene temper- 
ature change. We do not discuss that readout approach 
in this study. Instead, we treat graphene devices where 
the resistance (at the readout frequency) depends only 
weakly on temperature^^ so the temperature is read out 
by measuring the Johnson noise. Because dR/dT is very 
small, reading out the temperature change by measuring 



Johnson noise gives better sensitivity than readout using 
the resistance change with a finite bias current. 

The expected theoretical form^^ of Gcph = 4:AT,T^ 
is seen in the two experiments that measured electron- 
phonon coupling at low temperatures^'^'^^ (above 2 K). 
S measures the strength of electron-phonon coupling. 
Values of S with moderate electron density n w 2 x 
10^^ cm~^ are reported to be E = 70 mW/K^'m^ for 
graphene on Si02^^ and, with n w 10^^ cm~^, E ~ 

0. 5-2 mW/K'^m^ for graphene on boron nitride (BN);^'^ 
the values of E on BN arc smaller by a factor of up to 
140, so the choice of substrate material may affect E. 
E scales with electron density as ^/n for electron densi- 
ties far from the charge neutrality point (as in the case 
of these measurements); this does not explain the differ- 
ence in measured values of E. Effects of electron elastic 
scattering may significantly reduce the electron-phonon 
coupling in such graphene samples below 1 K.^^ Since 
there are no measurements of Gcph below 2 K, we take 
the measured form to compute Goph- 

A. Near-Equilibrium Device Noise 

We first consider the noise in the graphene device at 
or near equilibrium, where T Tq. There are two dom- 
inant sources of noise. The first is due to intrinsic en- 
ergy fluctuations in the device. These fluctuations can 
be thought of as energy carriers (i.e., photons, phonons) 
entering and exiting the graphene device at random."^" 
The second source of noise is the accuracy limit with 
which one can measure the temperature of a device us- 
ing Johnson-noise thermometry.'^^ 

1. Intrinsic Noise 

The rms intrinsic energy fluctuation of the detector, 
SEintr, is due to the intrinsic thermodynamic fluctua- 
tions. For measurements using the thermal response 
bandwidth corresponding to the thermal time constant, 
T = G/G, this is given by'^^ 

SEi^tr = {kBT'Cy/^ (2) 

STintr = (kBTyCy/^. (3) 

Here C is the heat capacity at temperature T; at Tq, G = 
Go and G = Go . 5Tintr is the rms equivalent temperature 
fluctuation. We find in the following section that readout 
noise is larger than this intrinsic noise. For that case, the 
optimum averaging time is indeed Tavg = G/G. 

The intrinsic energy resolving power (often called the 
energy resolution) i?intr is usually defined using the en- 
ergy fluctuation (full width at half maximum, fwhm), 
which is equal to 2v21n2 • SEinti ~ 2.35 • ^i^intr, where 
^Eintr is the rms fluctuation. Thus, at T = Tq- 

i?intr(fwhm) « 0.42 • E/^kBTgCo. (4) 
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However, for the measurement scheme considered, we 
read out the temperature of the device, not the energy. 
While 5T and 6E give equivalent information in the lin- 
ear range, we find that sensitive detection of THz photons 
occurs far beyond near-equilibrium range. As a result, 
the resolution in Eq. (4) is the relevant figure of merit 
only for operation in the linear range. A criterion based 
on temperature measurement will be developed in Sec- 
tion II B that is relevant for operation beyond the linear 
range. 

2. Amplifier Noise 

We now address the issue of amplifier noise and focus 
specifically on readout of the temperature using the emit- 
ted Johnson noise. This readout measures the emitted 
noise power over a time Tavg- The increase in the average 
emitted noise power at frequencies / < k^T /his given by 
P{t) = kB^T{t)B, where B is the observing bandwidth. 
We assume that the noise emission is characterized sim- 
ply with an instantaneous electron temperature, ^^^"^^''^^ 
that the detector resistance is impedance matched to the 
amplifier (achievable with n w 10^^ cm~^ and a wide 
graphene flake), and that we can use the low- frequency 
limit of the Johnson noise. '^'^ Thus, the signal measured 
in the absence of noise is given by 

S = [ " kBAT{t)Bdt, 
Jo 

where B is the output bandwidth and AT{t) = T{t)-TQ. 
If we take Tavg to be equal to the thermal time constant, 
T, this becomes: 

S = ksAT^.^BT, (5) 

where ATavg = i /o"(T(t) - To)dt. 

The rms accuracy with which the electron temperature 
T can be measured via Johnson noise in a time r is'^^ 

Ta + T 

^Trcadout = , , (6a) 

V Bt 

where Ta is the noise temperature of the amplifier. "^^ Near 
equilibrium, T « Tq and Eq. (6a) is equivalent to an rms 
energy fluctuation of 

(^i^rcadout = Cq ^. -— . (6b) 

V Bt 

The optimum values for B and Ta require careful 
consideration, and are constrained by amplifier choice. 
For the Johnson-noise readout, large bandwidth is desir- 
able if phonon cooling, Goph, dominates. This reduces 
(^Tlcadout- However, because Gphoton is proportional to 
B and r = C/G, once Gphoton > Geph, further increas- 
ing the bandwidth has limited benefit. To provide spe- 
cific examples of amplifiers, we will consider two ampli- 
fiers for the remaining calculations. The first (referred to 
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TABLE I. Summary of amplifier characteristics for the para- 
metric and SQUID amplifiers. 



as amplifier A) is a hypothetical near-quantum-limited 
parametric amplifier (paramp) which has a center fre- 
quency /o, bandwidth B, and noise temperature Ta that 
have all been scaled down by a factor of approximately 
10 from published data. "^^ Although the paramp specifi- 
cations listed in Table I have not been demonstrated, the 
measured results at frequencies approximately 10 times 
larger provide encouragement for their possible realiza- 
tion. The second (amplifier B) is a SQUID (Supercon- 
ducting Quantum Interference Device) amplifier with its 
measured performance."^^ The amplifier parameters are 
summarized in Table I. 

The total resolving power is limited by both intrinsic 
thermodynamic fluctuations and the accuracy limits of 
the temperature measurement: 

SEtot = Go ■ STtot 
Taking a large enough heat capacity to remain in the 
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FIG. 2. Prediction of energy resolution for a graphene detec- 
tor operating in the linear range with E = 0.5 mW/K''m^; for 
a 1-THz photon, this is for Co > 7 x 10"^° J/K, which corre- 
sponds to Of < 0.1 for To = 0.1 K. The readout is performed 
using amplifier A. The inset shows a normalized histogram 
of recorded signals from single-photon detection events and 
from sampling the baseline, with _Rtot = 1. This corresponds 
to a = 0.14. 

linear range, which we define as AT < 0.1 To, we find 
that the resulting calculated fwhm energy resolution, 

Rtot (fwhm) = 0.42 • E/6Etot, (7) 

is less than 1 with amplifier A. In Fig. 2, we plot this 
energy resolution in the linear range and also present the 
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normalized histograms that would result from counting 
1-THz photons with i?tot = 1 f^nd from sampling the 
baseline with no photon events. From this figure we see 
that it is not possible to clearly distinguish between the 
ensemble of zero photons, centered on signal, and the 
one photon histogram centered on 1. For the realistic 
case where the number of zero-photon events is much 
greater than the number of onc-photon events, the zero- 
photon peak would be much larger and the overlap of the 
two histograms would be worse. We therefore conclude 
that we must consider non-linear operation. Such opera- 
tion increases AT more than it increases i5Treadout, which 
remains the dominant source of broadening. 



B. Non-Equilibrium Device Noise 

In the previous section, we considered the device noise 
with no incident photons. However, there do not ap- 
pear to be any set of device parameters for T > 0.1 K 
that provide good sensitivity and keep operation in the 
linear range. We thus need to consider operation be- 
yond the linear range, for which the arrival of a pho- 
ton will significantly heat the graphenc. The temper- 
ature rise is computed from the electron energy, given 
by U — (7/2) AT^ at low temperature,'^'' where T is 
the electron temperature. Thus, for photon absorption, 
AU = E = i^/2)A{T^^^^~T^). 
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FIG. 3. Results from calculation of T(t) at Co = 2x 10"^^ J/K 
and B = 150 MHz with E = 0.5 mW/K^m^ after absorbing a 
1-THz photon. Also shown is the fit of the temperature to an 
exponential function, with To = 0.1 K and AT{t) calculated 
as discussed in the text. The effective time constant extracted 
is Toff = 460 ns. 

To determine 5i?intr outside the linear range the de- 
tailed time evolution of the temperature T(t) is con- 
sidered. In the linear range T{t) is given by AT{t) = 
Tit) - To - iE/Co)expi-To/t) for t > 0; tq = Cq/Gq. 
For large AT, the time decay is not a simple exponential 
function. The initial time decay at Tpeak is fast if phonon 
emission is dominant (r oc for phonon cooling), but 
the later decay of AT{t) back to Tq occurs more slowly. 
Additionally, for large output bandwidths, much of the 



cooling is done through emitted Johnson noise. In order 
to determine the electron temperature as a function of 
time, T{t), (shown in Fig. 3) we solve numerically 



?-7AT^ 
dt dt 

= -P{T) 

- -I]A(T4 - T^) - k^B{T - To) 
-iATdT 



-dt = 



(8) 



The computed T{t) — Tq is then fit to an exponential 
decay to determine an approximate effective time con- 
stant. Toff. Tavg is determined from this solution as 
shown in Fig. 3 using the weaker electron-phonon cou- 
pling, S = 0.5 mW/K^m^. Using Eq. (3) we can directly 
calculate the intrinsic thermal fiuctuations. The temper- 
ature T used in Eq. (3) is given with reasonable accuracy 
by (Tavg +2o)/2, as the temperature fiuctuations are due 
to phonons or photons leaving the graphene (at Tavg) or 
entering from the substrate (at Tq). The heat capacity 
of the graphene is given by C(Tavg). We use the aver- 
age electron temperature Tavg to specify T in Eq. (6a). 
The prediction of the temperature fiuctuations for this 
average temperature is shown in Fig. 4 as a function of 
Co, for To = 0.1 K with E = 0.5 mW/K^m^. In Fig. 4 
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FIG. 4. Prediction of rms temperature fiuctuations when 
excited by a 1-THz photon for To = 0.1 K and E = 
0.5 mW/K^m^. The readout calculation uses amplifier A, 
described in the text. Tavg is the average temperature of the 
graphene over TeS following the arrival of a 1-THz photon. 
For Co = 2 X 10"^^ J/K, Tavg - To = 0.48 K. 



we see that (^Treadout > i^Tntr- Thus, even though our 
evaluation of (5Tintr may not be exact, the readout noise 
dominates and determines <5Ttot, as seen in Fig. 4. The 
rms fluctuation widths in Fig. 4 apply for absorption of 
individual 1-THz photons. 

In Fig. 5, we plot the equivalent temperature fluctua- 
tions for the near-equilibrium case; this is relevant for the 
photon detector when there are no absorbed photons. In 
both Figs. 4 and 5, Tes depends on Cq. The Toff used to 
calculate each value of (5Ti-oadout and Tavg is different for 
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each heat capacity and is determined from fits analogous 
to the one shown in Fig. 3. The same values are used in 
Figs. 4 and 5. 
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FIG. 5. Prediction of rms near-equilibrium temperature 
fluctuations witli no photons for To — O.I K and E — 
0.5 mW/K''m^. The readout is performed using amplifier 
A. 

The energy resolving power will be evident if one plots 
a histogram of photon absorption events for a large en- 
semble of events. The histogram records the detected 
temperature rise for each absorbed photon. We denote 
this detected temperature rise as ATdct = 'S'/(fcBi?Tavg), 
with the measured value of S for each count. The en- 
semble average of the T^ct values is Tavg and the fwhm 
width of the distribution is 2.35 • STtot since STtot as 
defined above is an rms value. We plot this normal- 
ized distribution of single-photon counts in Fig. 6 for 
A^dot = 2dct ^ Tq. The functional form is 



Counts (X exp 



(^dct -^avg)'' 



^ "-'tot 



(9) 



For these plots, we choose Cq = 2 x I0~^^ J/K. This 
gives a larger ATavg than is obtained for larger values 
of Co. We do not consider a smaller heat capacity as 
it would lead to more significant heating of the electron 
system, potentially allowing the high energy tail of the 
hot electron distribution to diffuse over the energy gap of 
the superconducting contacts. With a carrier density of 
10^^ cm~^ this value of Cq = 2 x 10~^^ J/K corresponds 
to an area approximately equal to 4.5 /im^. Somewhat 
larger values of Cq (and A) would have similar perfor- 
mance. 

We also plot in Fig. 6 the predicted distribution if 
we sample the detector output with no input photons, 
using the same amplifier and Teff as employed for pho- 
ton detection. These are the "zero-photon" events, cen- 
tered around ATjot = 0. In the figure, the curves are 
normalized to the same height. However, in practice 
there will be many more zero-photon than one-photon 
events in the applications for which these devices are be- 
ing considered. ^'^^ Assuming one samples at a rate equal 
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FIG. 6. Normalized histograms with Co = 2 x 10"^^ J/K for 
photon counts using (a) amplifier A and (b) amplifier B. (c) 
shows the predicted photon count distribution if the larger 
value for E is assumed, with amplifier A. The dashed line 
in each figure indicates the estimated threshold temperature, 
ATdot.min, that is necessary to meet the single-photon detec- 
tion requirements described in Section III. 



and Tavg ~ 0.5 /zs, we expect there to be approximately 
200 more zero-photon events than one-photon events. 
For these plots, we define the temperature resolution as 
-Rtot.m (fwhm) = Aravg/(2.35 • STtot), where m can be 
or I to denote the zero-photon or one-photon case; i5Ttot 
is the rms value. 

Figure 6 clearly demonstrates that resolution is re- 
duced upon the arrival of an incident photon. However, 
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even the one-photon resolution is greater than is shown 
in the inset of Fig. 2. In that figure, the zero-photon and 
the one-photon histograms have the same width. With 
a small Co, we find that the resolution is improved for 
single-photon detection and is significantly improved for 
the zero-photon histogram. 



III. APPLICATION OF GRAPHENE-BASED 
DETECTORS 

The previous section makes clear that the tempera- 
ture resolution of a graphene calorimeter is not signal 
independent. After a photon has been absorbed, the fluc- 
tuations are significantly larger. The increase in fluctu- 
ation likely precludes the use of a graphene calorimeter 
for high-resolution THz spectroscopy or other methods 
of photon detection which require fine energy resolution. 
However, for single-photon counting, the outlook is pos- 
itive. 

One can use graphene as a single-photon "click" de- 
tector in which a threshold temperature rise ATdot.min is 
set. This threshold would be set above the main distri- 
bution of the equilibrium (zero-photon) fluctuations. We 
set ATdot.min to exclude almost all events with a smaller 
signal. This reduces the dark count rate below the speci- 
fied maximum, but may also exclude some actual (single- 
photon) events. The threshold needs to take into account 
the much larger number of zero-photon events than is dis- 
played in the normalized histograms of Fig. 6. 

To obtain a negligible dark count rate with the para- 
metric amplifier in Fig. 6a (amplifier A), one would need 
to set a threshold value for ATdot,min of 0.16 K. This 
threshold is indicated by the dashed line in Fig. 6a 
and it demonstrates that, with this choice of ampli- 
fier, nearly all single-photon events will be counted and 
almost all of the zero-photon events will be excluded. 
With the SQUID amplifier in Fig. 6b, a similar restric- 
tion on dark counts would require a higher threshold 
(ATdct.min = 0.36 K) and some single-photon events 
would be excluded. Worse performance is found by as- 
suming the stronger electron-phonon coupling (shown in 
Fig. 6c). A threshold ATdot,min = 0.4 K would need to 
be set to prevent excessive dark counts, excluding nearly 
half of the single-photon events. 

We conclude that a weaker value for electron-phonon 
coupling is necessary for an efficient detector using the 
Johnson noise readout. However, if the resistance change 
due to photon absorption is large, as would be the case 
with superconducting tunnel contacts, the predicted per- 
formance can be better. In this case the intrinsic energy 
fiuctuations could dominate and one might approach the 
intrinsic energy resolution (fwhm) shown in Fig. 2. 



IV. CONCLUSIONS 

Graphene detectors that use a Johnson-noise readout 
of detector temperature can count single THz photons 
with reasonable energy resolution. The design must care- 
fully take account of the effects of photon heating on 
the performance, and should employ graphene with weak 
electron-phonon coupling. It will be necessary to use 
temperatures below 1 K and low-noise amplifiers with 
nearly quantum-limited sensitivity. These requirements 
are largely necessitated by the Johnson-noise readout 
and the fast thermal response time of graphene. With 
a temperature-independent resistance (such as that of 
graphene), measurement of the electron temperature is 
difficult. Thus, the implementation of such a detector 
will require real care. 

Recent experiments^"^'^^ have shown that for some 
graphene samples the electron-phonon thermal conduc- 
tance can differ by up to two orders of magnitude. Mech- 
anisms that can decrease the electron-phonon emission 
rate in graphene are also under study.^^ Reduced phonon 
emission would increase the sensitivity for both single- 
photon and power detection. Thus, direct experimental 
tests and a more thorough understanding of how elec- 
trons interact with the lattice in real graphene samples 
will directly benefit the development of ultrasensitive 
photon detectors, and is one of the important current 
research challenges. 
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